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Fano-Kondo effect through two-level system in triple quantum dots 
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1, Komukai Toshiba-cho, Saiwai-ku, Kawasaki 212-8582, Japan 

We theoretically study the Fano-Kondo effect in a triple quantum dot (QD) system where two 
QDs constitute a two-level system and the other QD works in a detector with electrodes. We found 
that the Fano dip is clearly modulated by strongly coupled QDs in a two-level system and a slow 
detector with no interacting QD. This setup suggests a new method of reading out qubit states. 



I. INTRODUCTION 

Quantum dot (QD) systems have been providing op- 
portunities to probe a wide variety of many-body ef- 
fects of electronic transport properties in microelectronic 
structures. The Kondo effect, one of the main charac- 
teristics of these systems, is a result of quantum corre- 
lation between localized spin in QD and free electrons 
in electrodesi. Recently, the Fano effect, which appears 
as a result of quantum interference between a discrete 
single energy level and a major electronic system, has 
also attracted the interests of many microelectronics re- 
searches. Typically, a dip structure can be observed in 
conductance plotted as a function of an energy level of a 
side QlM^. 

A T-shaped QD system consists of two QDs, in which 
the first QD is connected to a source and a drain and 
the second QD is set beside the first QD. This system is 
well suited to investigation of both the Fano effect and 
Kondo effect (Fano-Kondo effect) and has been theoreti- 
cally analyzed by many authors^i^ii^. There are a variety 
of parameter settings in a T-shaped QD system. Wu et 
al. assume that there is an infinite on-site Coulomb in- 
teraction in the first QD (hereafter we call it 'detector 
QD'), and no interaction in the second side-QD^. On 
the other hand, Giiclii et al. assumes an infinite on-site 
Coulomb interaction in the second side-QD and no inter- 
action in the detcctor-QD^. In the former case, Kondo 
peak is dominant, whereas in the latter case, Fano dip 
structure is dominant in the transport properties of such 
a system. Tanaka et al. assume infinite on-site Coulomb 
interactions in both the detector-QD and the side-QD, 
showing that Fano dip structure appears in the middle 
of Kondo resonance peak^. 

Here we theoretically investigate the Fano-Kondo ef- 
fect in a triple QD system, depicted in Fig. 1. When 
two QDs are coupled, the energy levels mix, and bond- 
ing and anti-bonding states are formed. Thus, it is ex- 
pected that the current that flows through the detector- 
QD reflects these electronic states when coupling ratio 
iclt'd is sufficiently large, and the Fano dip would thus 
be modulated due to this two-level system. We can also 
regard this setup as an extra impurity connected to a 
T-shaped QD system, or a charge qubit closely attached 
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FIG. 1: Schematic plot of QD system. QDs a and b constitute 
a two-level system that is coupled to QD d only which is con- 
nected to the electrodes. We consider two cases with differing 
on the strength of on-site Coulomb interaction in QD d. In 
case I, there is no on-site Coulomb interaction in QD d. In 
case II, there is infinite on-site Coulomb interaction in QD d. 



to a QD detector. Coupled QDs can be used as a charge 
qubit^ dOd^i^^d^ . Charge qubits are usually capacitively 
coupled to detectors that are switched on only after quan- 
tum computation is finished. Results of quantum com- 
putation are readout through the field-effect (Coulomb 
interaction) depending on the positions of electrons in 
the charge qubit^ d^d^d^d'^ . If transfer of an electron be- 
tween a charge qubit and a detector is not prohibited 
(td 7^ 0) as in the present setup, detector current is ex- 
pected to be changed by the Fano-Kondo effect. This is 
an alternative method for reading out the quantum state 
of the charge qubit, although this method might change 
the quantum state of qubit more than the previous one, 
because detector electrons go through the charge qubit 
directly. To operate a coupled QD system as a charge 
qubit, we have to maintain the electronic state of the 
two QDs near the degeneracy point of charging energy, 
or hold the total number of excess electrons in the charge 
qubit to one. These requirements lead to stronger con- 
straints on this system. Here, as a first step, we consider 
only an effect of a two-level system without discussing 
qubit behaviors. 

A. W. Rushforth et al. experimentally investigated the 
Fano-Kondo effect when there are two side QDs near the 
electronic transport regiortii. They showed that Fano 
dip is modulated by the path of electrons through two 
side QDs. Compared with a simple T-shaped QD sys- 
tem, there are more parameters in this system due to the 
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existence of extra QD. Correspondingly, the Fano-Kondo 
effect in this system is expected to show more variety of 
dip behaviors than that in the T-shaped QDs as shown in 
RefJ^ii^. To make the problem transparent without loss 
of generality, we assume that there is a single energy level 
in each QD and that the two energy levels of the QD a 
and QD b coincide and correspond to gate voltages that 
arc applied to those QDs. For convenience, hereafter we 
call QD a and b a "two-level system" . We mainly investi- 
gate the change of conductance as a function of the gate 
bias. 

We use slave boson mean field theory (SBMFT) with 
the help of nonequilibrium Keldysh Green's functions. 
The formulation of SBMFT is very useful and a good 
starting point to study the transport properties of a 
strongly correlated QD system. Thus, this method is 
widely used, although it is only usable at a lower tem- 
perature (T) region than the Kondo temperature Tr -^^'^^ , 
that is T < Tk. Here, we investigate two cases differing 
in the amount of Coulomb interaction in QD d. In case I, 



there is no Coulomb interaction {Ud = 0) and in case II, 
there is an infinite Coulomb interaction (U^ = oo). The 
case I experimentally corresponds to a large QD d and 
the case II corresponds to a small QD d. It is expected 
that we can see clear Fano effect in case I, and a Kondo 
resonance peak strongly changes Fano dip in case II. Ac- 
tually, it will be shown that, in order to see the effects of 
bonding- antibonding states in the detector current, cou- 
pling between the two side QDs should be strong and the 
detector should be sufficiently slow only in case I. 



II. FORMULATION 

The Hamiltonian for the case II is different from that 
for the case I in that the former has an additional con- 
straint on QD d. The mean field Hamiltonian for the case 
II is described in terms of slave bosons bai {ai = a, b, d) 
as: 



Ci—L,Rka,s cti—a.b.ds ai—a.b.d \ s / 

+ ^ T.if^-Mfbs + flhbifas) + ^ Y^^flMlhs + flM\fds) + E ^ll^^i.sb\hs + flMu^s) (1) 
s s a=L.R k^,s 
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where Ek^ is the energy level for source (a = L) and 
drain (a = R) electrodes. Ea, Eb and Ed are energy 
levels for the triple QDs, respectively, tc, td and Vq, are 
the tunneling coupling strength between QD a and QD 
b, that between QD b and QD d, and that between QD 
d and electrodes, respectively, c^^s and /q^^ ai'c annihi- 
lation operators of the electrodes, and of the triple QDs 
(ai = a,b,d), respectively, s is spin degree of freedom 
with spin degeneracy N; here we apply N = 2. A„j 
is a Lagrange multiplier. We take = and 
Ea-^ = Ea-^ + Xai as mean field parameters. The Hamil- 
tonian for case I is similar to H^^^^ except that = 
and 6<i = 1 in Eq.([l|). 



By applying a conventional procedur e i to this 
Hamiltonian, we derive Green's functions for QDs, for 
example, G'^^{uj) = Dab/Boo etc., where Dab = - 
Ea){io - Eb) - il, Boo ^ DabBr - {lu - Ea)\id\^ and 
Br=Lu-Ed + iT with ic = tcbabl/N and id = tdbdbl/N 
(sec Appendix). Here, F^ = 2TTpa{EF)\Va\'^ is the tun- 
neling rate between a electrode and QD d with a density 
of states (DOS), pa{Ep), for each electrode at Fermi en- 
ergy Ep. F = (Tl + r_R)/2 and we assume F^ = Tb,- 



Four self-consistent equations for case I are given. 

icY.(fLfas)+Xa\ba\''=0, (2) 

s 

icJ2ifaJb^) + idY.{fLfbs)+h\bb\' = 0, (3) 

s s 

J2(fLfo..^) + KA'^l. («i = a,6), (4) 

s 

For case II, we have to add two more equations regarding 
constraints on QD d with Va — Vabd/^/N- 

i*dY.^flks)+Y. Y.v:{ci^jds)+Mbd\'=Q (5) 

E</i/<i.> + Ibdl' = 1- (6) 

s 

To see the Fano-Kondo effect, we calculate conduc- 
tance G = diL/dV at zero bias V = 0. Source 
current II is calculated from Green's functions as 
lL = fJ ^^^FflF4/i(c.) - fniu;)) with /4^) ^ 
[exp{{uj + eV)/T) + l]-^ and /^(tj) = [exp(a;/T) + 
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Then, we obtain a conductance formula: 



G 



Bn 



(7) 



The ratio tc /td compares the internal coupling strength 
in a two-level system with that between the two-level 
system and the detector. Here, we regard the case where 
ic/td = 5 as a strongly coupled two- level system and 
the case where tc/td = 1 as a weakly coupled two-level 
system. Because TL^n/iVL +^r) (= F ; we assume 
= Ffl) can be approximately estimated as a tunneling 
rate through the detector, if F/i^ is sufficiently large, the 
electron that flows through QD d is so fast that it cannot 
detect the oscillation of an electron in the coupled QDs 
a and b. In contrast, if T /td is small, the electron that 
flows through QD d can observe the oscillation between 
bonding and antibonding states. Therefore, we call a 
detector with large F/t^ = 2 a fast detector, and one 
with smaller T /td = 0.4 a slow detector. Tk is estimated 
as Tk ~ j:)(^-^\Ei-Ep\/T ^ ^ Q^^ p jg ^ bandwidth), 
when we assume that D = 20td, \Ed\ < OAtd, F > OAtd 
and Ep ~ 0. 



III. NUMERICAL CALCULATIONS 
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FIG. 2: Conductance G as a function of an energy level of the 
two- level system Ea(= Eb) with a strong coupling {tc/td — 5) 
in case I {Ud = 0). (a) T/td = 0.02 and (b) T/td = 0.2 for a 
fast detector (T/td = 2). (c) T/td = 0.02 and {d)T/td = 0.2 
for a slow detector {T/td = 0.4). Comparison of (a) with (c) 
shows that the dip structure is modulated in a slow detector 
when energy level Ed is changed. Ef = 0. Numerical results 
are obtained in 10~* relative error range. 



Here, we numerically investigate the transport proper- 
ties of our triple QD system for case I {Ud = 0) and case 
H {Ud^^). 



A. Case I: Ud ^ 

Figure [2] shows a conductance G for strongly coupled 
two- level system {tc/td — 5) in case I. For a fast detector 
with T/td = 2 (Fig. [H^a) and (b)), clear Fano dip can be 
seen in both low temperature {T/td — 0.02) and high 
temperature {T/td = 0.2). However, we cannot see the 
effect of the change caused by the bonding-antibonding 
state, even if we change Ed- When F is reduced as shown 
in Figl2] (c), we can see that the Fano dip is changed 
from a simple dip form to an asymmetric form at low 
temperature. The modulation of the structure in the 
Ed <Q region is considered to be an effect of the bonding- 
antibonding oscillation in the two-level system. These 
behaviors can be understood from DOS, which is derived 
from pd{'^) = -ImG^^(a))/7r : 



t^^{LU - Eaf - 2tj{uj - Ed){uJ - Ea)Dab 



Instead of a single dip in the numerator in a T-shaped 
QDs^, Dab in the numerator of Eq.® indicates dou- 
ble dips, that comes from bonding-antibonding state of 
side-QDs. Figure [3] shows four DOSs, which are dis- 
played corresponding to FigO We can see that there 




FIG. 3: Density of states Eq.((8ll for a strong coupling {tc/td = 
5) in case I {Ud = 0). (a) T/td = 0.02 and (b) T/td = 0.2 for 
a fast detector {T/td = 2). (c) T/td = 0.02 and {A)T/td = 0.2 
for a slow detector {T/td = 0.4). Ed = -Q.2td-EF = 0. 



Dlb\Br? 

(8) 



is a single peak for Ea/td — —2.4, closely distributed 
three peaks for Ea/td — —1.2, and separated three peaks 
for Ea/td = 1.2. By comparing Figs. [3l[a) and (c) with 
Figs.[3^b) and (d), respectively, we can see that temper- 
ature dependence on the shape of DOS is weak. Tem- 
perature slightly shifts the peak position. In contrast, F 
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FIG. 4: Mean field values as a function of Ea(= Ei,): solutions 
of the self-consistent equations in case I. (a) fast detector, (b) 
slow detector. T = 0.02td. = ~Q.2ti. U = 0.05. Ep = 0. 
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FIG. 5; Conductance G as a function of Ea{= Ei,) for a weak 
coupling {tc/td = 1) in case I [Ud = 0). (a) T/td = 0.02 and 
{h)T/td = 0.2 for a slow detector (F/td = 0.4). Figure (a) 
shows that the dip structure is modulated for a slow detector 
when energy level Ed is changed. Ef ~ 0. 



dependence is large, similar to a Breit-Wigner function. 
The DOS of the slow detector (Fig. [3] (c)(d)) is sharper 
than that of the fast detector (Fig. [3](a)(b)). F mainly 
determines the width and the maximum value similar to 
T-shaped QD system. Because Eq.([7]) is expressed as 
G = -i2e/h) J dLjpd{ij)zdriirL/iTH + TL)dfL{oj)/duj, 
conductance reflects an energy region ±T around Ep in 
DOS. From dfL{i^)/duj — > S{uj) when T 0, at low tem- 
perature, conductance can detect the change of closely 
distributed three peaks in the DOS if those peaks are 
sufficiently sharp such as a slow detector. On the other 
hand, the present results show that the peak structure 
for a fast detector is not sufficiently sharp to detect the 
three peaks. At higher temperature, conductance inte- 
grates peak structures in DOS, resulting in the single-dip 
structures shown in Figl^] (b) and (d) . These are reasons 
of higher sensitivity of the slow detector at low temper- 
ature. 

Figure |4] also shows the difference between a slow de- 
tector and a fast detector. Oscillations of mean field val- 
ues for a fast detector (FigH] (a) ) are weaker than those 
for a slow detector (Figg] (b)) at E'^ < 0. 

Figure [S] shows the case of weak coupling {tc/td = 1) 
in a slow detector. Fano dip modulation exists at lower 
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FIG. 6: Conductance G as a function of an energy level of the 
two- level system Ea{~ Et) for a strong coupling {tc/td ~ 5) 
in case II {Ud = oo). (a) T/td = 2 and {h)r /td = 0.4. Ed = 
Ef = 0. 



temperature, but is smaller than that of Figure [SJc). At 
the higher temperature region of T = 0.2td, we can ob- 
serve a small and calm dip structure for both fast and 
slow detectors. 



B. Case II: Ud = oo 

In this case, full self-consistent equations are applied 
so that the total number of electrons in QD d is less than 
one, and we have to solve the six self-consistent equa- 
tions from Eq.dll) to Eq.®. Without QD a and QD 
6, a Kondo peak is expected. In the case of T-shaped 
QDs, the Fano dip appears in the middle of Kondo peak 
conductance^. In the present case, the two energy levels 
in the side-QDs couple with the Kondo peak and pro- 
duces the complicated conductance structure shown in 
FiglH Here, overall structure of G does not change even 
in the region oi Ed < 0. 



IV. DISCUSSION 

As noted in the introduction, the present model can 
be applied to several situations. If we look at QD a as 
an extra trap site for the T-shaped QD system composed 
of QD b and QD d with both electrodes, the present cal- 
culation shows that a simple dip structure is most likely 
to be observed, and when the extra trap site is strongly 
coupled to the side QD b and the energy level of QD d is 
sufficiently low, the modulation of Fano dip could be ob- 
served. Experimentally, this would happen because it is 
not easy to exclude all trap sites from the QD system due 
to the difficulty of nanofabrication^. The present calcula- 
tions also indicate that when there are many additional 
trap sites, the Fano dip structure would be modulated 
further. 

To construct a charge qubit^ii^ in the present case, 
the total number of electrons in QD a and QD b should 
be one. However, in the present setup, we only apply 



5 



less than one electron in each QD. Thus, we need more 
restriction to hold one electron through the QD a and QD 
6, for example, by introducing an exchange interaction 
such as that discussed in Reii^. 

Although SBMFT is believed to be valid in previous 
works for one or two QDs, it has not yet been proved 
that this approximation is valid for many QD systems 
in the context of the Fano-Kondo effect. The effect of 
fluctuations around mean fields should be investigated in 
the near future. 



V. CONCLUSION 

We have studied the transport properties of a triple QD 
system, in particular, emphasizing the Fano-Kondo ef- 
fect through two-level system. We have used slave-boson 
mean field theory to describe quantum interference be- 
tween electrons in discrete energy levels and free elec- 
trons in electrodes. We found that the Fano dip struc- 
ture is modulated from simple dip to clear asymmetric 
form, when the coupling strength between the QDs in 
two-level system is large, using a slow detector without 
on-site Coulomb interaction in the detector QD. We also 
found that the Fano dip is strongly changed by Kondo 
resonance peak when there is an infinite Coulomb inter- 



action in the detector QD. 
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APPENDIX A: GREEN'S FUNCTIONS 

Green's functions are obtained based on equations of 
motion for the time ordered Green's functions. Fourier- 
transferred diagonal Green's functions for QDs are given 
as GM = [{iJ - Eb)Br - \id\^]/Boo and Gl,{u;) - 
[{u! — Ea)Br]/ Bqo othcr than GJJ^(a;) expressed in the 
main text. We also have off-diagonal Green's func- 
tion by using relations (w — Ea)Gab{'^) = tcGbb{^) and 
{u;-Eb- \ia\y{cj - Ea))Gab{uj) = i^GM{uj)- GauA^) 
and Gk^k'^ {i^) {a = L, R) are obtained by {cj — E'^ — S;, — 
\id?/[u: ^Eb - Itrfl - Ea)]]Gdk^ {u) = Vo^liu - EkJ 
and {uj - Ek^)Gk^k'^{uj) = Sk^k'^ + VaGdk^i^^), with 
= Y.a J2k^ I^«IV('^ - ^kj- After applying analytic 
continuation rules, we obtain lesser Green's functions in 
the current formula in the main text. 
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